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In quantum mechanics, systems can be described in phase space in terms of the Wigner function 
and the star-product operation. Quantum characteristics, which appear in the Heisenberg picture 
as the Weyl's symbols of operators of canonical coordinates and momenta, can be used to solve the 
evolution equations for symbols of other operators acting in the Hilbert space. To any fixed order in 
the Planck's constant, many-body potential scattering problem simplifies to a statistical- mechanical 
problem of computing an ensemble of quantum characteristics and their derivatives with respect to 
the initial canonical coordinates and momenta. The reduction to a system of ordinary differential 
equations pertains rigorously at any fixed order in ft. We present semiclassical expansion of quantum 
characteristics for many-body scattering problem and provide tools for calculation of average values 
of time-dependent physical observables and cross sections. The method of quantum characteristics 
admits the consistent incorporation of specific quantum effects, such as non-locality and coherence 
in propagation of particles, into the semiclassical transport models. We formulate the principle 
of stationary action for quantum Hamilton's equations and give quantum-mechanical extensions of 
the Liouville theorem on the conservation of phase-space volume and the Poincare theorem on the 
conservation of 2p forms. The lowest order quantum corrections to the Kepler periodic orbits are 
constructed. These corrections show the resonance behavior. 

PACS numbers: 02.30.Hq, 02.30.Jr, 02.70.Ns, 05.30.-d, 05.60.Gg, 25.70.-z 



I. INTRODUCTION 



The deformation quantization [l|, 0, d, 0, IE El has been the focus of renewed interest last decades. It uses the Weyl's 
association rule [l], || to establish the one-to-one correspondence between phase-space functions and operators in the 
Hilbert space. The Wigner function Q appears as the Weyl's symbol of the density matrix. Refined formulation of the 
Weyl's association rule is given by Stratonovich (tJ. Groenewold [i| proposed the quantum-mechanical description of 
the evolution in phase space and introduced into the formalism the concept of star-product of phase-space functions. 
The quantum evolution is determined by the skew-symmetric part of the star-product, known as the sine bracket 
or the Moyal bracket [H, @. The Moyal bracket represents the quantum deformation of the Poisson bracket, being 
essentially unique Q. The phase-space quantum dynamics keeps many features of the classical Hamiltonian dynamics. 
The formulation of quantum mechanics in phase space and the star-product are reviewed in Refs. [y, ITol [Til . 
14, 15 1 . The Stratonovich version Q of the Weyl's quantization and dequantization is discussed in Refs. 
17 . 18. [T9L [20|. Wigner functions have found numerous applications in quantum many-body physics, kinetic 
21 , 22 I , collision theory, and quantum chemistry [23|, [24| ■ Transport models created initially for needs of 




quantum chemistry in order to describe chemical reactions were modified and extended for modelling heavy-ion 

collisions 0, H, M, E3, H, M, H SH H m . 

First attempts to exploit specific properties inherent to the *-calculus have recently been made towards calculation 
of higher order terms of the ^-expansion of the Bohr-Sommerfeld quantization rules [371 ] and determinants of operators 
with occur in quantum field theory at one loop [38l l39l . l4(i| . A diagrammatic method for calculation of symbols of 
operator functions in form of the series expansion in h has been proposed [U, \^ . 

Transport models in heavy-ion physics are designed for a phenomenological description of complicated reaction 
dynamics of nuclear collisions. There exist several well established transport models such as the Boltzmann-Uehling- 
Uhlenbeck (BUU) approach [HEl, (Relativistic) Quantum Molecular Dynamics (R)QMD [H, H3, [H, [H, M [3JJ 
or Anti-symmetrized Molecular Dynamics (AMD) [43|, [44j. These approaches can be derived from quantum field 
theory [25j, |45j, |46|, |47| in semiclassical limit and contain basic quantum features such as final state Pauli blocking 
for binary collisions of fermions. Numerical solutions are realized by propagating test-particles (BUU) or centroids 
of wave packets (QMD, AMD) along classical trajectories in phase space. In the case of AMD, wave packets are 
anti-symmeterized according to their parameters. Transport models provide a sound phenomenological basis for 
description of a variety of complicated nuclear phenomena. Quantum interference effects are, however, beyond the 
scope of these models. The problem of intrinsic consistency of approximations remains under the discussion and 
stimulates further developments (see e.g. [H, |H, |49[ and references therein). 

The most striking feature of the semiclassical transport models is obviously the phase-space trajectories along which 
the particles or wave packets of the particles are supposed to propagate. While scattering of n classical particles can 
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be processed by means of conventional computer programmes, the evolution of many-body wave functions represents 
a field-theoretic problem that can currently be approached neither analytically nor numerically. Any implementation 
of the consistent quantum dynamics should obviously keep trajectories as an attribute allowing to access many-body 
scattering problems numerically. 

The notion of phase-space trajectories arises naturally in the deformation quantization through the Weyl's transform 
of Heisenberg operators of canonical coordinates and momenta [l^, S3, [H3| . These trajectories obey the Hamilton's 
equations in the quantum form and play the role of quantum characteristics in terms of which the time-dependent 
symbols of operators are expressed. In classical limit, quantum characteristics turn to classical trajectories. The 
knowledge of the quantum phase flow is equivalent to the knowledge of the quantum dynamics. 

In this work, we report the semiclassical expansion of quantum characteristics for many-body scattering problem, 
based on the method [l^, [HI, 113, H3] j and provide tools for calculation of average values of time-dependent physical 
observables and scattering cross sections. 

We show that to any fixed order in h it is sufficient to work with quantum characteristics, provided equations 
of motion and rules for calculation of the time dependent average values of physical observables are modified. The 
quantum evolution problem becomes thereby entirely identical with a statistical-mechanical problem. We hope that, 
to a fixed order in h, the method of quantum characteristics captures basic quantum properties of many-body potential 
scattering being numerically effective the same time. 

We thus propose a self-consistent non-relativistic quantum mechanical approach for solving many-body potential 
scattering problem of spin-zero particles. The method is based on semiclassical expansion of quantum characteristics 
and star-functions of quantum characteristics in a power series over the Planck's constant. The evolution equations 
represent a coupled system of first-order ordinary differential equations (ODE) for quantum trajectories in the phase 
space and for associated Jacobi fields. 

The outline of the paper is as follows: In the next Sect., we overview the Weyl's association rule and the star-product 
using the method proposed by Stratonovich [7J. Sect. Ill is devoted to discussion of quantum characteristics. We 
investigate transformation properties of canonical variables and phase-space functions under unitary transformations 
in the Hilbert space. The role of quantum characteristics concurs with the role of classical characteristics in solving 
the classical Liouville equation. Quantum characteristics are physically distinct from trajectories of the de Broglie - 
Bohm theory [5l[ and trajectories appearing in the semiclassically concentrated states (3|| [36| . 

In Sect. Ill, we derive the principle of stationary action for quantum Hamilton's equations, extend to quantum 
systems the Liouville theorem on conservation of the phase-space volume by the classical phase flow, and propose the 
quantum counterpart of the Poincare theorem on conservation of 2p forms. 

Sect. IV is devoted to the semiclassical expansion of a star-function around an ordinary dot-function. Specific 
features of this technique applied to star-functions of quantum characteristics are discussed in Sect. V. We expand 
quantum characteristics in a power series in H and derive a system of coupled ODE for components of quantum 
characteristics of the h expansion and for the associated Jacobi fields. We construct in particular the Green function 
for the lowest order quantum correction to classical phase-space trajectory in terms of the Jacobi fields. 

The lowest order quantum corrections for the classical Kepler periodic orbits are constructed in Sect. VI. 

Numerical methods of solution of many-body scattering problems using the h expansion and calculation of average 
values of physical observables in the course of quantum evolution are discussed in Sect. VII. In Sect. VIII, we 
summarize results. 



II. WEYL'S ASSOCIATION RULE AND THE STAR-PRODUCT 



Classical systems with n degrees of freedom are described within the Hamiltonian framework by In canonical 
coordinates and momenta £ l = (g> , q n ,Px, ■ ■■,p n ) which satisfy the Poisson bracket relations 



(HI) 



with the matrix / 



-E n 

E n 



(II.2) 



where E n is the n x n identity matrix. In quantum mechanics, canonical variables are associated to operators of 
canonical coordinates and momenta f = (q 1 , q n ,pi, ...,pn) € Op(L 2 (M. n )) acting in the Hilbert space, which obey 
the commutation relations 



(II.3) 
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The Weyl's association rule extends the correspondence <-> f to arbitrary phase-space functions and operators. 

The set of operators f acting in the Hilbert space is closed under the multiplication of operators by c-numbers and 
summation of operators. Such a set constitutes the vector space V . The elements of its basis can be labelled by 
canonical variables £*. The commonly used Weyl's basis looks like 

®(0 = (2nhr5 2n (t - y) = / — exp(--ry fe (C - y) fe ). (II.4) 



The Weyl's association rule for a function /(£) and an operator f has the form [3] 

/(0 - Tr[B(Ofl, (H.5) 

In particular, = Tr[2$(£)y 1 ]. The function /(£) can be treated as the coordinate of f in the basis 23(£), Eq. pi.5[) as 
the scalar product of 33(£) and f. 

Alternative operator bases and their relations are discussed in Refs. 0,(3]. One can make, in particular, operator 
transformations on 2$(£) and c- number transformations on £\ 

The Weyl-symmetrized functions of operators of canonical variables have the representation [52j 

f=/(^^))' (H.7) 

where the subscripts indicate the order in which the operators act on the right. 

The set of operators is closed under the multiplication of operators. The vector space V is endowed thereby with 
an associative algebra structure. Given two functions /(£) = Tr[58(£)f] and g(£) = Tr[<B(£)g], one can construct a 
third function 

/(O*3(O = rr[B(0fe] (II.8) 
called star-product Q . It is given explicitly by 

f(0*9(0=f(0exp(^V)g(t), (II.9) 

where 



d d 



is the Poisson operator. 

The star-product splits into symmetric and skew-symmetric parts 

/*5 = /°.9+y/A<?. (11.11) 

The skew-symmetric part / A g is known under the name of Moyal bracket. 

The Weyl's symbol of a symmetrized product y n y 12 ...f s simplifies to a pointwise product of canonical variables 

rr[Q3(0p ( V' 2 ...? ia) ] =C 1 -C", (H.12) 

which is explicitly symmetric with respect to permutations of the indices. The symmetrized product of In Hermitian 
operators u l is associated to the symmetrized star-product of real phase-space functions = Tr[2$(£)u 1 ]: 

rr[B(Ou (<a u <a ...u i ' ) ] = m (h (£) o u l2 (£,) o ... o it ls) (0- (H.13) 

The o-product is not associative. The order in which it acts in Eq. pi.13p is, however, not important since the indices 
are symmetrized. 

The Weyl's association rule can be reformulated in terms of the Taylor expansion. Consider the Taylor expansion 
of a function /(£): 



/(o-EasFSsM*-^- (,U4 » 
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For any function /(£) one may associate an operator fa- 



riJIL,, x is (1I15) 



The simple calculation 



-f 1 ...^ (IL16) 



l d s f(o) 

s=0 

^ 1 S /(O) 
s—Q ^ 

shows that /t(£) = /(C) an d fr = f- The Taylor expansion over the symmetrized products of operators of canonical 
coordinates and momenta gives the association rule equivalent to Eqs. (|II.5[) and (|II.6[) . So, one can write f = /(y). 

The average values of a physical observable described by an operator f is calculated as trace of the operator product 
fr where r = r + is the density matrix, Tr[r] = 1 or, equivalently, by averaging the phase-space function /(£) over the 
Wigncr function 

W(0=Tr[«8(0t]. (11.17) 

The Wigner functions is normalized to unity 

d 2n £, 



(2?rS) 



-W(0 = 1. (11.18) 



If f <~> /(£) and r <-> W(£), then 



Tr ^ = / (0/(0 * no = I ^mw(o. (n.19) 

The star-product can be replaced with the pointwise product 0, EH, • 

Not every normalized function in phase space can be interpreted as the Wigner function. The eigenvalues of density 
matrices are positive, so given r, one can find a Hermitian matrix z 1 / 2 such that r = ti^t^. For any W(£) there 
exists therefore W 1/2 (g) such that W(£) = W 1/2 {Q * W 1/2 (£). 

For a pure state r = \tp >< the Wigner function becomes W(£) =< ip\%$(£)\ip >. Its value is restricted by 
—2™ < W(£) < 2™ provided \tfj > has a finite norm (56| . 



III. QUANTUM CHARACTERISTICS 



This section is devoted to studying the transformation properties of the Weyl's symbols of operators under the action 
of unitary evolution operators. One-parameter set of unitary transformations applied to the Heisenberg operators 
of canonical variables generates, under the Weyl's association rule, phase-space trajectories. The knowledge of these 
trajectories is equivalent to the knowledge of quantum dynamics. In particular, time-dependent symbols of operators 
are functionals of such trajectories. In this sense, the phase-space trajectories play the role of characteristics. The 
interest to the quantum characteristics is connected with the special role of trajectories in the transport models. 

In Euclidean space, we can construct basis {ej} and study two kind of transformations: Passive ones where are 
transformed and other vectors do not transform, and active ones where are not transformed and other vectors are 
transformed. These two views are equivalent. In what follows, active transformations are considered, the operator 
basis 25(£) remains fixed. This corresponds to the Heisenberg picture where evolution applies to operators associated 
to physical observables. 



A. Unitary transformations under the Weyl's association rule 



Consider unitary transformation of an operator f — * f = it + fil where it + il = = 1. The operators of canonical 
variables are transformed as 

f* -f f = ii+fii, 
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whereas their symbols as — > = Tr[58(^)il + p l iX]. Define 



= Tr [®(C)il + p l it]. (III.l) 



The associated transformation of function /(£) has the form 

HO - AO = Tr[«8(0fl - Tr[Q5(0it + fil] 



7i 



s! d^K.-d?* 



s=0 

= /MO)- (in-2) 

The ^-product can be replaced with the o-product. The o-product is not associative. The order in which it acts is, 
however, not important due to symmetrization over the indices. As a consequence, semiclassical expansion of /(*£*(£)) 
around f(u(£)) involves even powers of H only. In general, f(*u(£)) — /(oit(£)) ^ /(«(£)), while /(*«(£)) = /(«(£)) 
provided is a linear function in £\ 

The antisymmetrized products y^ 1 ...y l2s l of even numbers of operators of canonical coordinates and momenta rep- 
resent c-numbers. They are left invariant by unitary transformations: 

ii+r^---y i2s] ii = f [ll ---f J2sl - (ni.3) 

In phase space, equations (|III.3|) look like 



2 ) (2s)! 



Y^(-yp^..,p*°-M> (IH.4) 



where summation runs over permutations of indices («i«2---i2s) and (— ) cr = ±1 depending as the ordered set (iii2---*2s) 
constitutes even or odd permutation of (1, 2, 2s). In particular, 

«'(£) A = f A f = -i^ ■ (HI.5) 

One may associate to real functions u 1 ^) Hermitian operators f = u l (y). If functions M l (£) obey Eqs. (|III.5jl . 
operators f l obey the commutation rules for operators of canonical coordinates and momenta 

[u i (t),ui(t)] = [f,?] = -ihiv. (III.6) 

There exists therefore a unitary operator il which relates f and y\ [68| Applying Eq. (|III.2[) to the product fg of two 
operators, we obtain a function /(C)*<7(C)Ic=*«(£,t) associated to the operator il + (fg)il and function /(*w(£))*<?(*u(£)) 
associated to the operator (il + fil)(ll + gil). These two operators coincide, so their symbols coincide also: 

/(C) *g(Qk=Mt) = f(MO)*a(MO)- (ni.7) 

The star-products in the left- and right-hand sides act on £ and £, respectively. 

Equation (|III.7j) shows that one can compute the star-product in the initial coordinate system and change the 
variables £ — > C — or equivalently, change the variables £ — > C — * u (£) an d compute the star-product. Equation 

pil.7[) applies separately to the symmetric and anisymmetric parts of the star-product. 

Equation (|III.7[) makes it possible to calculate the star-product in new unitary equivalent coordinate systems. 
The functional form of equations constructed with the use of the summation and the star-multiplication operations 
remains unchanged in all unitary equivalent coordinate systems. The star-product is not invariant under canonical 
transformations in general [l9l |. 



6 



B. Quantum phase flows generated by one-parameter sets of unitary transformations 



A one-parameter set of unitary transformations in the Hilbert space is parameterized by 

H = e X p(-lfrr), (III.8) 
n 

with S) = $j + being the Hamiltonian. The functions defined by 1|) acquire a dependence on the parameter r, 
so one can write r). They specify quantum phase flow which represents quantum analogue of the classical phase 
flow [57j] • In virtue of Eq. (|III.2[) , the evolution of symbols of Heisenberg operators is entirely determined by ti*(£, r). 

We keep the conventional term 'canonical transformations' for functions preserving the Poisson bracket. Phase-space 
transformations which preserve the Moyal bracket (|III.5|) are referred to as 'unitary transformations'. This is consistent 
with the fact of discussing the continuous unitary operators (jIII.8l) . Continuous groups of unitary transformations 
represent the non-trivial quantum deformation of continuous groups of canonical transformations. As we distinguish 
between the Poisson and Moyal brackets, we have to distinguish between the canonical and unitary transformations. 

The relationship between the canonical transformations, which are not elements of continuous groups, and trans- 
formations in the Hilbert space appears to be more involved [HI, |j54| . 

The energy conservation in the course of evolution along quantum characteristics implies 

H(0 = H(Mtr)) (III.9) 

where 

H(Q = 7V[B(Oi5] (IIL10) 

is the Hamiltonian function. Quantum characteristics can be found from the Hamilton's equations in quantum form 
(see e.g. [H) 

^u i ^,r) = {C,H(0}\ c=Mtr) (ni.ll) 

with initial conditions 

u l (£,0)=f. (111.12) 

The substitution £ = *tt(£, r) to F' 1 ^) — H(()} leads to a deformation of the classical phase flow. If the ★-symbol 
would be missing, the quantum effects could be missing also. 

The phase-space velocity <9m 1 (£, t)/<9t depends on u*(£,t) like in classical mechanics and on derivatives of m 1 (£,t) 
with respect to f , as a specific manifestation of quantum non-locality. 

Functions which stand for physical observables evolve according to equation 

/&r) = Tr[<B(£)iI+fit] - /(*«(£, r),0), (111.13) 
whereas the Wigner function does not evolve W(£,t) = W(£, 0). Functions satisfy equation 

-^m, r) = /(£, t) A H(0 = /(C, 0) A H(0k= MM (111.14) 

which is the Weyl's transform of the quantum-mechanical equation of motion for Heisenberg operators. 
The formal solutions to Eqs. pil.lip and (|III.14[) in the form of series expansions in r are given by 

00 s 

u ^< r ) = E^T (4eAffK))Aff(0)A-^)). (HI.15) 

S ^ s 

OO g 

/(*«(£, r)) = ^ (-(( f(0 A H(0) A ff (0) A ...H(0). (111.16) 

s=0 ' 

s 

The transformation of canonical variables to order O(r) has the canonical form, since A H (^) = {£% H(£)}. The 
second order in r gives deviations from the canonical transformation. The infinitesimal transformations generate 
canonical or unitary finite transformations according as the composition law is t±), T2) = u 1 (£,ti + T2) or 

n), T2) = it l (£,Ti + T2), respectively. The infinitesimal transformations of functions which are not linear in 

are not canonical and cannot result from a c-number transformation of the coordinate system. 

If an operator 21 commutes with Sj, its symbols is conserved in the sense of A(t;) — A(-ku(£,T)), similarly to 
Eq. (|III.9|) . In the Schrodinger picture, the Wigner function of a stationary state commutes with S) and obeys Ws (£,) = 
W a(-k u($ —t)). For a harmonic oscillator, m 1 (^,t) depend linearly on and coincide with the classical trajectories 
0,111, [59[. In such a case, the ★-symbol can be dropped to give Ws(£.) — Ws(u(£, — t)). 
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C. Green function in phase space in terms of characteristics 



We wish to establish a connection of the phase-space Green function [6CJ, |61| with quantum characteristics and to 
derive Eqs. pil.lip using the Green function method. 
Combining Eq. (|II.5[) and Eq. (|II.6[) . one gets 

f = / i0r m)Trm)f] - (IIL17) 

In Eq. piI.13|) , we substitute on place of f the right-hand side of Eq. 17|) and obtain 

M, r)= J C, r)/(C, 0) (HI-IS) 

where /(£,0) = /(0 and 

£>(& C, r) = D(C, £, -t) = Tr[mm +f B(Cn = Tr[B(C)il«8(0il+], (HL19) 
with D(£, C, t) being the Green function. The use of Eqs. pil.17p and piI.19[) for 11 = U 2 iXi gives 

J ^^ D ^'M D ^^) = D ^ 1 '^ Tl - + ^)- ( IIL2 °) 

The Green function can be treated as the Weyl's symbol of il + Q3(C)H in the basis 58(£) or the Weyl's symbol of 
il*8(£)il + in the basis 25(C). T ne orthogonality condition 

Tr[®(0®(0] = {2TrK) n 5 2n {i - C) (111.21) 
follows from Eqs. pl.5|) and (|II.6[) . Using Eq. (jIII.2|l . we obtain 

D(i, C, r) = (2TrK) n 6 2n (*u($, r) - C) = (2wh) n ^ - *u(t, ~r)). (HI.22) 

The evolution equation for the Green function can be found by calculating the time derivative of Eq. 19|) and 
applying the Weyl's transform to the commutators of f) and il +s S(()ii or il<B(£)il + : 

—D{t C, r) = D(£, £ t) A H(0 = -£>(£, £ r) A (C). (111.23) 

One can replace in Eq. piI.18[) /(£ 0) by f* and /(£, r) by u l (£, r), take the time derivative and apply Eq. pil.23[) . 
We get 

j^D& C, r) * (C A 2T(0) = y 7^7^ C, r) (C A # (0) . (IH.24) 
Using the identity Q AH(Q = {£, H(()} and the explicit form of the Green function (|III.22p we arrive at Eqs. pil.lip . 

D. Principle of stationary action for quantum Hamilton's equations 

The quantum Hamilton's equations can be derived from the principle of the stationary action. The action S is the 
integral 

^ lr dT l (4* iU<( *' T) ** ( * ,T) ~ F( ^' T) ' (m ' 25) 

We compare trajectories u(£,t) coming out £ at t\. The action is stationary provided 2n functions u z (^,t) obey 
equations pil.lip . 
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The stationary action principle can be reformulated in more familiar terms for Heisenberg operators of the canonical 
coordinates and momenta: 



S= drTr 

' Tl 



We recall that H(i) = Sj is the Hamiltonian. Using 

k i ,H(r)] = ih{?,H(t)}\i=f, 

one arrives at the Heisenberg equations of motion 



(111.26) 



(111.27) 



(111.28) 



and, further, at Eq. (IIII.ll[) upon the Weyl's transform. 



E. Quantum counterpart of the Liouville theorem 



The Liouville theorem states conservation of the phase-space volume by the classical phase flow. The Poincare 
theorem of Hamiltonian dynamics suggests conservation of the 2p-forms under the canonical transformations (see e.g. 
[57J). We describe quantum analogues of these remarkable theorems. 




FIG. 1: Quantum phase flow applied classically to region Q. in phase space. The phase-space volume is not conserved with 
time. The region Q. at r = is specified by its index function in n (£)- The region Q at r ^ is specified by the evolved index 
function «n n («(C, — t)). 

Let n be a region in phase space. Its index function is defined by 

®® = {o, tin. ( IIL29 ) 

The function Iq™(£) can be associated to a region ft smaller than it is allowed by the uncertainty principle. It does 
not fulfill the condition of divisibility satisfied by the Wigner functions: For every W(£) one can find Wi/ 2 (£,) such 
that W(£) = Wi/2(t) *W / i/2(C)- The index function does not belong to the set of continuous functions, C°°(T*M. n ), 
whose derivatives all are continuous. There are no reasons to restrict with the set C°°(T*R n ), however, since infinite 
sequences of continuous functions have limits which are not continuous functions in general. Applying the Weyl's 
association rule to in n (£)i we obtain a Hermitian operator i^". 

The phase-space volume occupied by f2 at r = is calculated classically 

S2 "(°) = / TTT^^iO = Tr[i2fl. (III-30) 



Quantum phase flow preserves the Moyal bracket rather than the Poisson bracket and does thereby not transform 
canonically. The evolution of Q along quantum characteristics, respectively, does not preserve the phase-space volume: 

r d 2n f 



K n (t) = I J^^na, -r)) + E^(0). (111.31) 
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The variance appears at the order 0(h 2 ). 



^(★tV(t-T)) 




FIG. 2: A visual demonstration of quantum evolution of desity distribution associated to index function iff (£) in phase space 
at t = 0. The boundary of region Q diffuses in the course of quantum evolution. The integral of the evolved index function 
in n (*u(f, — t)) over phase space is conserved. 

One can make a positive statement, however. Implicitly, we treated the transform classically substituting if"(£,) — > 
io — t)). Quantum particles, however, do not move along quantum characteristics, as must be clear from the 
composition law 

n + r 2 ) = u (*u(& n), r 2 ). (ni-32) 

The boundary of f2 does, naturally, not move along quantum characteristics also. It experiences quantum fluctuations. 
In the spirit of quantum mechanics, one can speak of conservation of average volume only. The quantum evolution is 
formally expressed through the star-product transform £ — + *it(£, — r), so we have to write 

-r)) = Tr[ili^il+]. (HI.33) 



Equation 



— Tr[fl = (III34) 



is valid for any Heisenberg operator f G Op(L 2 (]R")) including i^ 1 . 

Remarkably, it incorporates the phase-space volume conservation within the framework of the Groenewold-Moyal 
dynamics: 

S 2 "(r) = £ 2n (0). (111.35) 

In the limit of h — » 0, r) — > C (£, T ) where c(£,t) are solutions of the classical Hamilton's equations. One has 
therefore limfi_>o £ 2 ™(r) = lim^o S 2rl (r), and we recover the classical Liouville theorem. Equation (|III.34[> and its 
consequence E)q. (|III.35|) comprise the quantum- mechanical analogue of the classical Liouville theorem. 

F. Quantum counterpart of the Poincare theorem on conservation of 2p-forms 

The quantum analogue of the Poincare theorem on conservation of the 2p-forms under canonical transformations 
follows from Eq. (|III.34|) also. 

Let H,2p be a 2p dimensional submanifold in phase space. It can be specified by constraint equations G a (0 = for 
a = 1, 2m where m = n — p. 

The volume element of fi 2 p is given by [13] 

J2nt 2m 

= 7^( det ^m G b (0}) 1/2 (^r II S ^ C &)- (IIL36) 

^ * C— 1 
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The index function Iq 1 ^) of f^p can be taken to be 

2m 

vf (0 = (aet{g a (o,G b (m) 1/2 (^rT[s(G c (0)4 n (o- (ni.37) 

0=1 

where in (£) is an index function in phase space non-singular across the constraint submanifold. The delta- functions 
force i^ p (^) to vanish outside of £l 2p - 

The function (|III.37j) can be associated to an index operator i^f by the Weyl's association rule Eq. (|II.6|) . The 2p 
dimensional volume S 2p (r) has the form of Eq. pil.33|) with jQ n (£) replaced by ijf(£). 

The quantum phase flow, in virtue of Eq. (IIII.34j) . preserves the 2p volume 

E 2p (r) = £ 2p (0). (111.38) 

For p — n, we recover Eq. pil.35l) . 

In classical mechanics, 2p- forms take values on vectors d^2, d&p to give oriented volumes of 2p-dimensional 
parallelepipeds. In the course of evolution, infinitesimal parallelepipeds remain infinitesimal parallelepipeds. 

In quantum mechanics, we can construct infinitesimal parallelepipeds at r — and calculate their volumes also. As 
distinct from the classical case, however, quantum evolution does not keep index functions i 2 |'(^,0) within the set of 
index functions. For t ^ 0, «q P (£,t) does not correspond to any region in phase space, although phase-space integral 
of Iq (£, t) is equal to volume (or area) of the initial parallelepiped at r = 0. There is only one instant of time r = 
when the formalism shows the clear geometric sense. 



IV. SEMICLASSICAL EXPANSION OF /(*«(£, r)) AROUND /(«(£, r)) 

The evolution problem for quantum systems can be split in two parts: First, we look for quantum characteristics 
which are solutions of Eq. (|III.llj) and, secondly, use equation /(£, r) — /(*u(£,t),0) to calculate time-dependent 
symbols of operators. The key problem is to have an efficient algorithm to calculate functions of £ = *u(£, r). The 
★-arguments appear in the quantum Hamilton's equations (jlll.lip and the evolution equation for functions 141) . 

We give semiclassical expansion of /(*m(£, t)) around /(m(£,t)). Function /(£) can be represented as the Fourier 
transform 

^=j^M^ k e)m. (iv.i) 

It is sufficient to learn how to calculate exp(*[7) where 



U = l -l lk u\i 1 T). 



Using Eq. pTL2|) . we obtain [6£ 



where 



exp(*C7) = (1 + h 2 c 2 + h 4 c 4 + 0(h 6 )) exp(J7) (IV.2) 



c 2 = -^-{2UUP 2 U + 3UP 2 U), (IV.3) 
48 

c 4 = ^^(90(^P 2 ^)^P 2 ^ + 60(^P 2 ^)P 2 ^ + 48(^l/P 2 C/)[/P 2 ?7 + 45(C/P 2 C/)(C/P 2 C/) 

+ m{UUP 2 U){UP 2 U) + 20(UUP 2 U){UUP 2 U) + 30{UUP 2 U)P 2 U) 

+ T-^( 6UUUUpi U + teUUUP^U + 30([/£/P 4 C/)i + 40(PPP 4 £/) 2 + lbUP A U). 
11520 

Owing to combinatorial factors, P acts like V: In particular, we have UP 2 U = lfP 2 U. In the expression UUP 2 U, 
V 2 acts one time to each U on the left and two times on the right. More sophisticated expression (UP 2 U)P 2 U is 
computed by acting P 2 two times on (UP 2 U) and two times on U . Expressions of the type ABP 2 C are computed 
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by acting V one time on A, one time on B and two times on C, e.g., (UP 2 U)UP 2 U stands for A — (UP 2 U) and 
B = C = U . The formal description uses notations 



A(0' h - 



d s A(Q 
(-yi i ^...I i ^A(0, jl ... js . 



(IV.4) 
(IV.5) 



We have 



AP 2 B 


= A, kl B< kl , 


ABP 2 C 


= A tk B tl C> k \ 


AP i B 


= A.jjir" 1 . 


(ABP 4 C) 1 


= A,JS u C' jU . 


(ABP 4 C) 2 


= A ti B >jkl C' ijkl , 


ABCP i D 


= A.H.CmI)-' 11 . 


ABCDP 4 E 


= A.n.Cjj),/-:"' 1 



(IV.6) 



The diagram technique [4l| can be useful for high-order ^.-expansions. Our calculation pV.3p is in agreement with 

m. 

The lowest order quantum correction to f(*u(£, r)) can be found using Eq. (|IV.2p . We replace r\i — > —ihd/du l and 
U — > u l (£)d/du l and obtain 



/M«>t)) = M,t)) 



T) tl v>(Z, T), m u k (t ry m f(u(Z, r)) >ijk 
^u < «,r), w ^(^T)> w /(««,r)),« + 0(fi 4 ). 



(IV.7) 



The next-order correction can be found using Eqs. (|IV.2p and pV.3[) . The derivatives of f(u(£, r)) are calculated with 
respect to u l . 

As an application, we consider the semiclassical expansion of the Weyl's symbol of the projection operator = 
S(E — Sj). The result is in agreement with Ref. [3]. The series expansion pV.2[) with U = —H(£)/T can be used for 
calculation of the Weyl's symbol of the finite-temperature density matrix. 



V. SEMICLASSICAL EXPANSION OF QUANTUM CHARACTERISTICS 

In this section, we discuss general properties of the quantum characteristics, specific for the semiclassical expansion. 

A. Semiclassical expansion of quantum Hamilton's equations 



The first step in solving evolution problem is to construct quantum characteristics using Ecis. (|III.ll[) . We make an 
expansion in a power series of the Planck's constant: 



Ut) = X> 2 X(£,t). 



(V.8) 



s=0 



Here, u l (£, t) is the classical trajectory which starts at r = at a point £\ The right-hand side of quantum Hamilton's 
equations 1 1|) . F l (() = {£', -ff (C)}> i s a function of £ = t), so we have to use a power series expansion similar 
to (HVJl) : 



F i (M^T)) = J2h 2s Ft(u (^T),...,u s ^,r)). 



(V.9) 
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Given u(£, r), functions FJ(uo(£,t), u s (£, r)) are known. i^ s l depend on derivatives of t) for < r < s with 
respect to In particular, 

F > ) = ^(uo), 

fj(uo,ui) = uj(£,r)F>ob ~ ^ 3 oK^),™4K,t),„4K,t)^"P( U o) jH 

- ^ugCCTj.tatig^r)'""^^)^, (V.10) 

and so on. To any fixed order in the Planck's constant, quantum characteristics can be found by solving finite-order 
partial differential equations (PDE) 

Jx = j?|(u , ...,«,), cv.il) 

with initial conditions 

«o&o) = e, 

<(C,0) = 0, s> 1. (V.12) 

Given t) for < r < s, equations for t) simplify to a system of first-order ordinary differential equations 
(ODE). Such a circumstance allows to approach the problem recursively using numerically efficient ODE integrators. 

Classical trajectories Mq(^,t) appear at O(H ). To order H 2 , the first and second order derivatives of Mq(£,t) with 
respect to canonical variables are involved. We have therefore to consider not only propagation of points in phase 
space, like in classical mechanics, but also propagation of gradients 

at < r < s which affect high-order quantum corrections r) and enter semiclassical expansion of r)). 
We call such gradients generalized Jacobi fields (tensors) or simply Jacobi fields, since first-order derivatives, which 
determine stability of classical trajectories with respect to small perturbations of parameters, are known in classical 
mechanics as the Jacobi fields. We use the symplectic form T lk to shift the indices up and down according to the rule 
(|IV.5|) . e.g., x % — —P k Xk, x-i — —IikX k where Iui — ~I kl , so that Xiy 1 = —x l yi. In what follows, we discuss r = 
Jacobi fields, so the lower index r will be suppressed. 

B. Semiclassical expansion of energy conservation condition 

The energy conservation in quantum form (|III.9[) gives a sequence of conserved quantities 

H s (u (t t), ...,«.(£, r)) = |m 1 ^ ffMC, r)). (V.14) 
For s > 1 they vanish at r = 0, so that 

H s (u o ,...,u s ) = (V.15) 

provided H(£) does not depend on h. The term s = in Eq. (|V.14|) gives the energy conservation law for classical 
systems (cf. ([HO]) ) 

H($ = H(u (Z,t)). (V.16) 

The functions H s and F* depend linearly on u\ for s > 1: H s (uq, ... ,u s ) — u l s H{u )s + H s (uo, ...,it s _i) and 
F*(uo, u a ) — —r lk u J s H(uo)^j -\-Fj(uo, u s -i). Using equations of motion, one can remove from time derivative 
of H s , so that 

d 

—H s (u , ...,u s ) = G s (u , ...,u s _i). (V.17) 
or 
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In case of s = 1, the quantity 

1 



GiM = -^{J l i(^r)J^,r)J klm ^,T),H(O}H(u ), l3k 
1 



16 -{Jk(£. ^ & r), fl"(0}H(«o),« (V.18) 

remains constant G\(uo) = on classical trajectories. It does not depend on quantum corrections u*(£,t), although 
involves derivatives of Uq(£,t) with respect to £\ We have not found for Gi(uq) known counterpart in the classical 
Hamiltonian theory. [70| 



C. Green function for first-order quantum corrections in terms of Jacobi fields. 

Let us establish some useful identities for Jacobi fields. The energy conservation implies 

Taking derivative from this identity, one gets 

~d^ = — a* — ( v - 19 ) 

These are the six integrals of motion involving both the classical trajectory and its rank-two Jacobi fields. 
By writing the Hamilton equations in the form 

r)n i (P t) 

— = {u^,r),H(u (tr))} = {ufo,r),H(0} 

one gets 

^^ = Ji«,r)(-2«)^. (V.20) 

From canonicity of the classical hamiltonian flow one has 

J^T)I rs J k s ^,T) = -{ul^,r) > u^ ) r)}=I ik . (V.21) 
Equation (|V.21|) implies that —J r k (£,,T) is the inverse matrix of J* - (£, t), and so 

jr(Z,T)J£(t,T) = SI (V.22) 

J%foT)J?.(t,T) = SI (V.23) 

Multiplying Eq. (|V.20[) by J.- t (^,t) and performing the summation over i one gets Eq. (|V.19[) . These equations are 
therefore not independent. 

The evolution equation for the first-order quantum correction has the form 



du 



= -i ik d2H d %x )) u[ ^ T) + riM ^ (v - 24) 

where ^(uoi^r)) is the inhomogeneous part of Eq. (|V.10|) . For periodic orbits, the system of equations for u\(£, r) 
looks like a system of equations for coupled oscillators whose frequencies are time-dependent. The functions J l (u(£,, r)) 
play the role of an external time-dependent force. 
We look for solutions of equation (|V.24|) in the form 

4(£,t)=C s (£,t)J^,t). (V.25) 
Taking equations of motion (7.1) for the rank-two Jacobi fields into account, one gets 

J! s &t)-^C s (£,t) = J> (£,r)). (V.26) 
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Using Eq. (|V.23[) and the initial conditions u\(0,£) = 0, we obtain 

t*i(£,T) = r) f T J r s (r',0J r (Mr',0)dr'. (V.27) 

Jo 

The function 

G r s (r,r',0 = ~J r m (tr)J s m (r\O0(r-r') 
entering the right-hand side of Eq. (|V.27|) satisfies equation 

(*sf + "' d2H a7 t{ )r ) «> = *'< T - T '> < v ' 28 > 

and can be recognized as the Green function for the first-order quantum correction to the classical trajectory. 



VI. QUANTUM CORRECTIONS TO KEPLER ORBITS 

The classical Kepler problem is described in many textbooks. This problem has the exact solution in quantum 
mechanics also (see e.g. [63]). We construct the lowest order quantum corrections to Kepler orbits. 

In order to calculate Jacobi helds, one has to keep the initial canonical variables as free parameters. It means that 
we have to analyze Kepler orbits in arbitrary coordinate system K' . When all the derivatives over £ are calculated, 
one can pass to a coordinate system K where the orbits belong to a [x, y) plane and where the expressions simplify. 



A. Spherical basis 

In the Cartesian coordinate system, the basis vectors are defined by 

= (1,0,0), 
e y = (0,1,0), 
e z = (0,0,1), 

The Kepler problem is usually treated in the spherical basis. The spherical basis vectors parametrized by radius r, 
polar angle and azimuthal angle ip have the form 

e r = (sin^costpjSin^sin^jCOsfl), (VI. 29) 

e e = (cos#cos<y9, costfsin^, — sm9) — — - , (VI. 30) 

o9 

1 de r 

= (-sin^,cos<^,0) = -— — -. (VIM) 
sm9 dtp 

They are orthonormal and obey the scalar and vector product rules: 
for i,j = r,9,p. 

In the spherical basis, the particle velocity can be decomposed as follows: 

dre r . ■ de r . de r 

x = ^ = rer + r9 W + rip ^ 

= re r + r9eg + r sin 9e v (p. 



The Lagrangian of the system: 



™(f 2 +r 2 9 2 +r 2 sin 2 9p 2 ) + -. 
2 r 
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The canonical momenta are as follows 

p r = mr, (VI.32) 

p e = mr 2 9, (VI.33) 

p v = mr 2 tp sin 2 9. (VI.34) 

In the spherical basis, the canonical momenta have the form 

(p r sin 9 cos tp + — cos 9 cos p sin tp \ 

p r sin sin tp + ^ cos 9 sin 93 + cos tp | (VI. 35) 
p r cos 9 — ^ sin 9 



r rsm( 



The Hamiltonian function equals 



^ _ _Pr_ , , & (vi 36) 

2m 2mr 2 2mr 2 sin 2 6 r' 



The orbital momentum L = x x p can be represented as follows 



Pep 

The absolute value of the orbital momentum equals 



L=p g e v -^-e e . (VI.37) 
sm9 



7-2 2, p y 



sin 2 (9 

It is useful to parametrize pg and p v as follows: 



pe = Lsinx, (VI.38) 
p v = Lcosxsin9. (VI. 39) 



B. Euler rotation 



We wish to pass over to a coordinate system K' in which the particle trajectory belongs to the (x' , y') plane, while 
the orbital momentum is directed along the z'-axis. Furthermore, we require the perigelium of ellipse along which the 
particle moves be on the positive side of the x' axis. 

In the initial coordinate system, K, the Cartesian components of the orbital momentum can be found to be 

— sin x sin tp — cos \ cos $ cos <Pi 
sin x cos f — cos X cos 9 sin tp, 
cos x sin 9. 

The vector L determines plane (x',y') orthogonal to L. This plane intersects the plane (x,y) of the coordinate 
system K. The line of the intersection is parallel to e z x L. It forms an angle, <pE, with the x-axis. The first Euler 
rotation brings parallel to e z x L. In order to find <pE, we notice that cosipE ^x(g z x L) = (e^ x e z )L = — L y , 
sin(j)E a e y( e z x L) = (e y x e z )L = L x and (e 2 x L)i = e^kLk = (—L y ,L x ,0), so that (e 2 x L) 2 = L 2 — L 2 . The 
second Euler rotation brings e z parallel to L. The first pair of the Euler angles (^e^e) can therefore be found from 
equations 

— sin x cos tp + cos x cos 9 sin tp 

y/T— cos 2 x sin 2 9 

— sin x sin tp — cos \ cos 9 cos tp 

\/\ — cos 2 x sin 2 
cos x sin 9, 

\/l — cos 2 % sin 2 0. 



COS <f>E = 

sm(f)E = 

COS^e = 

sin Ae = 
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Suppose the third Euler rotation by the angle ipg brings the x-axis parallel to the major half-axis of the ellipse 
in the direction towards the perigelium. Let us compare the unit vector n = e r in the two coordinate systems. The 
Cartesian components of n are related as follows 

sin 9 cos <p 
\E\\ [ sin sin 
cos 9 

where \\E\\ is the Euler rotation matrix: 

cosV>£ sink's \ / 1 \ / cos(pE sin^s 

\E\\ — | — sinV'B cost/'b cos^e sin^^ — svacpE coscpE 
1 / \ -sintf B cosd E J \ 1 

and <fp is the phase set off, i.e., the angular distance from the perigelium at start of the motion (r = 0). We obtain 

cosip'p = sin 9 cos ip cos ipE cos 4>e — sin 9 cos ip sin ipE cos $e sin 4>e 

+ sin 9 sin ip cos ip e sin 4>e + sin 9 sin ip sinipE cos'&e cos 4>e 

+ shxtpE sin e cos 6>, 
sin^p = — sin 9 cos i/3 sin ipE cos — sin 9 cos <p cos tpE cos sin ^ 

— sin sin y> sin ipE sin 0g + sin 9 sin <p cos ipE cos i9e cos 0e 

+ cos ipE sin ??b cos 0, 
= sint^s sin6>sin(0£ — ip) + cos-Oe cos 9. 

The last equation is the identity. The first two ones give 

cos ip' = Pe sin 6 cos ipE H = cosflsin^, (VI. 40) 

sinyj' = — Pe sin9 sin t/>^ H = cos 9 cos tpE- (VI. 41) 



These equations can be solved for ipE to give 



cos^sin^ - sin x sine cos ^ , lrT 

cos^b = = E, (VI.42) 

V 1 — cos 2 x sin 9 

. cos6»cos<^' +sinxsin6>sin^' , 

smVB = = -. (VI.43) 

V 1 — cos 2 x sin 



The canonical momenta in the coordinate system K' can be found to be 

p r cos ip' p - j sin p' p 
p r sin ip' p + - cos p' p I , 


so we conclude (cf. IVI.35|) p' r — p r , p' e — 0, and p' = L. Of course, 



\E\\ 













>-( 















Ly 






L z 1 







C. Orbits in the coordinate system K 



We use the atomic units m = a = % = 1. The classical orbits in the coordinate space are known to be the ellipses. 
The value p (do mot mix with momentum!), eccentricity e, major half-axis a of the ellipse have the form (see e.g. 
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my- 

P = L\ 

e = + 2EL 2 = cos x- 



P 

a = 



sin 2 x 



Given these parameters are known as functions of the initial canonical variables £ l = (r, 9, ip,p r ,pQ,p v ), one can find 
the set off angle <p' from equation 

(VI.44) 



1 + cos x cos If' 



The Euler angle ipE is then fixed by Eqs. (|VI.42[) - (|VI.43|) . The Euler angles (4>e, #e, ^e) are then determined. 
In the coordinate system K' , the classical orbit has the form 

TV = r' r = — V - 7 , (VI.45) 

1 + cos x cos (p' T 

r T = r' T = a(l — cos x cos v) , (VI. 46) 

r = p3/2 ^-«>BXBin^ 
sin x 

, COSV — COS X .„S 

cos^ = VI.48 

1 — cos x cos v 

, sin x sin v 

smip > = * VI.49 

1 — cos x cos ^ 

Prr = Prr = & , (VI.50) 



Per = 0, (VI.51) 
p' VT = VP- (VI-52) 

The motion starts at v v which corresponds to r p Eq. (|VI.47j) and (p' p Eq. (|VI.44|) . The initial conditions are as follows: 

7T 

2 : 



the values ip' p and L = y/p are functions of £. 

Now, we are in a position to find the angular coordinates of the trajectory in the K system: 



sin 9 T cos ip T \ / cos 

sin^sin^,- = sin<^ 
cos 6> r / \ 

where 1 1 £7 1 1 1 is the inverse Euler matrix. The radii in K and K' coincide. The canonical momenta in K may be 
found with the help of Eqs.(jVL32) - fVPl) . 

We thus get the phase-space trajectories in the initial coordinate system: 

w o(£' T ) = ( r T, T , (pT,PrT,P9T,P<pr)- 

D. Jacobi fields in parametric representation 

The phase-space trajectories represented in parametric forms bring specific features in the calculation of the Jacobi 
fields. In our case, the orbits are parametrized by the parameter v which is a function of time r and initial canonical 
variables £*. The derivatives of the functions u % (y, £) which we have constructed explicitly involve the derivatives of 
v = v(€,t). 

Let t p (£) be time needed to come to £ from the perigelium (x' — r, y' — 0). We have 

r ^) =p 3/2 ^-cosXsin^ (VI53) 
sin x 
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The dependence on £ comes through the parameters p, \i ancl v v The parameter v p — v p (P) is determined from 
Ecis. rfvmi) . (TVL481) and (TVL491). 

At t > t p , we shift in Eq. (|VI.47j) the time variable r — > r + r p (£) and obtain 



, U \ 3/2" ~ cosxsmv _ 
sin % 



(VI. 54) 



According to the convention, the motion starts at r = 0. Now, the parameter r is independent of £*. Taking the first 
and second derivatives of Eq. (|VI.54[) . we obtain 



dr p (Z) = dTjvX) dTjvX) dv 

dp dp dv dp' 

d 2 r p {Q _ d 2 T(v,p d 2 T{v,j) dv d 2 T(v,p dv 

dpdp* dpdv d£ k dp*dv dp 

d 2 T(v,p) dv dv <9T(v,£) d 2 v 



dpdp° 



dv 2 d£ k dp 



dv dpdp 



(VI.55) 



(VI.56) 



These equations allow to find the derivatives J^r and 



d 2 U 

WW 



at fixed t. The dependence on r enters implicitly through 



Given the time dependence in the parametric form, one finds the total derivatives and the Jacobi fields in the 
parametric representation, accordingly: 



du(v,p) du(v,p) dv 

dp + dv dp' 
d 2 u l (v,p d 2 u l (v,Q dv d 2 u l (v,p dv 
dpdp* + dpdv d£ k + d£ k dv dp 
d 2 u l (v,p dv dv du L (v,p d 2 v 

dv 2 dp d^ + dv dpdp* ' 



(VI.57) 



(VI.58) 



The formulae simplify considerably if we set K' = K after taking the derivatives. The explicit parametric form 
of Ji l (pT) and J l ik(P T ) can be found using MAPLE. We give results for J l i(pr) for circular orbits only (r) = tt/2) 
where one can set v p = ip p = 0: 



\fi(t>T)\ 



cos v 








p 3 / 2 sin v 





2P 1 ' 2 (1 - cosv) 





cosv 








— 1/2 • 

p 'smv 





— 2p _1 sin v 





1 


2p x l 2 {cosv — 


1) 


p~ 1/2 (4 sinv - 3 v) 
2p~ 1 sin v 










cosv 








—p i smv 








cosv 




















1 



(VI.59) 



The rank-two and rank-three Jacobi fields which we have constructed satisfy for arbitrary values p and \ the initial 
conditions J l ^(0,p — 5\ and J^(0,£) = 0, equations of motion (|VII.1|) . the conservation laws Eq. (|V.20[) . and the 
canonicity condition (|V.21[) . 



E. Quantum orbits and resonance 



In general case, the functions u{(p r) are expressed as one-dimensional integrals. Here, we report results for circular 
orbits, where significant simplifications occur. 
The external force is found to be 

J 3 = p~ 5/2 (-l + cosv)(5cos 2 v - 7 cosv + 1), 
J 4 = 3p~ 3 (-l + cosv)(2cos 2 v-2cosv- 1), 
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FIG. 3: The lowest order quantum Kepler orbit coresponding to the classical circular orbit r = p = 110 in the atomic units 
m = a = h = 1 for < v < 8n. The motion starts at x = r and y = 0. 



while J 1 = J 2 = J 5 = J 6 — 0. The components i = 1,2, 5, and 6 vanish for all orbits (0 < 77 < 0). Also, for all orbits 
J3 _ p-5/2 and j4 _ p -3_ 

The external force is periodic. The matrix entering the left-hand side of Eq. (|V.28j) for circular orbits has the form 



, d 2 ?iMg,T)) 

du r du k 












1 




















p 





2p- 5 /2 














p- 2 
















2 p -5/2 






































Inspecting this matrix, we find that frequency of the external force coincides with the eigenfrequency of effective 
oscillator. This condition is sufficient for the resonance. 

The Green function can be constructed using Eq. (|VI.59[) . The first-order corrections become 

3 7 15 

u\ = — cos 3 v + 4 cos 2 v H — cos v — 5 H vsmv, 

4 4 4 

3 -1 f ~ ■ 4 • 13 • 2 15 \ 
u\ = p — 10 sin v cos v H — smv-\ sini/cos v H i^cosi' — v , 

\ 3 6 2 / 

4 —3/2 / ^ 2 o • r» • 1^ 

tti = p ' — cos v sm ^ — 8 cos v sin ^ + 2 sin ^ H z^cosi/ 

\4 4 



while other components vanish: u\ — u\ = u\ = 0. The corrections u\, uf, and u\ vanish for < r\ < also. The 
observed dependence of u\, u\ , and u\ on p holds for arbitrary 77 and ip' p . 

The amplitude of the quantum corrections increases linearly with time, which is the feature inherent to the reso- 
nances. 

In Fig. [3l we plot the radius r T = r + r% versus the azimuthal angle (p T — ipo + ipi in the polar coordinate system 
for p = 110 during the four z/-periods (0 < v < 8tt), the motion starts at v = ipo = 0. The value of p is taken to be 
large in order to fulfill the semiclassical conditions where the first quantum correction is supposed to be small. The 
region of validity of the solutions is restricted to |i4(£,t)| <S |uq(£,t)|. We see from the plot, that for p — 110 one-two 
periods is the upper admissible limit to which the first-order quantum characteristics can be extrapolated. This is in 
the obvious contrast with the operator methods where the semiclassical approximation usually starts to work from 
very low quantum numbers. 
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VII. MANY-BODY POTENTIAL SCATTERING 

Suppose we have an n-body system which evolves quantum-mechanically due to an interaction through a potential. 
The initial-state wave function is assumed to be known, so we can construct the initial-state Wigner function. 

A. Calculation of Jacobi fields using bunches of characteristics 

The classical trajectories m (£,t) are assumed to be constructed numerically. In order to evaluate Jacobi fields 
(9*w (£, t)/<9£ 11 ...d£ H to order 0(h 2 ), one may consider a bunch of characteristics coming out from a neighborhood of 
£*. Using the simplest finite difference method, 3 2 ™ such characteristics should be calculated and stored. 

The order 0(h 4 ) involves Jacobi fields of u (£,t) with t = 1,2,3,4 and Jacobi fields of u\(^,t) with t = 1,2. 
Accordingly, one has to construct a bunch of at least 5 2 ™ characteristics u 1 ((,t) and at least 3 2 ™ characteristics 
u\((,t) where C's belong to a neighborhood of £\ 

The appearance of a bunch of characteristics is a manifestation of non-locality in quantum mechanics. Moving 
up with the Planck's constant expansion, one has to generate increasing bunches of characteristics. The number 
of characteristics increases exponentially. In the potential scattering of two particles, we would have to store ~ 
10 6 characteristics, whereas in the potential scattering of four particles, we would already have to store ~ 10 11 
characteristics to order ti 2 . 



B. Calculation of Jacobi fields using ODE for Jacobi fields 

It is possible to reduce the amount of data generated and stored by computer. The main problem is to evaluate 
Jacobi fields. 

Remarkably, Jacobi fields can be constructed by solving a system of ODE: 

^(uoXfcJ^K.T), 

F>oWi-(£,T)JI-(£,T) +F l (u )jJ l jk (C,r) (VTI.1) 
with initial conditions 

4k(t>°) = °- ( m2 ) 

Equations (jVII.ip follow from (|V.11[) for s = 0. Equations (jVII.ip are sufficient for constructing characteristics 
to order 0(h 2 ). To find the first correction «1(£,t), it is necessary to calculate (2n) 2 functions 9u (£, t)/8£ j and 
(2n) 2 (2n + l)/2 functions <9 2 u (£, t)/9£ j 9£ . One has to store, therefore, 2n(n + l)(2n + 1) functions describing 
characteristics and their derivatives instead of 2n x 3 2n ones describing a bunch of characteristics. 

The method of the previous subsection requires computer time of about r c ~ 2n(2s + l) 2 ™, whereas the present 
method gives r c ~ (2n) 2s+1 to order 0(h 2s ). The number of characteristics and Jacobi fields in the present method 
grows exponentially with s. It means that for a few-body system, a high order H calculation is more effective in terms 
of bunches of characteristics, whereas for a many-body scattering and a low s the method based on the propagation 
of Jacobi fields is preferable. 

To order 0(H 4 ), one has to construct derivatives of m (£,t) up to the fourth degree. We thus have to store 2n 
functions u (£, r) and 2n X} s =i( s + 2n— l)!/(s!(2n — 1)!) Jacobi fields. The corresponding equations can be obtained 
by taking derivatives of the classical Hamilton's equations with respect to initial values of the canonical variables. 
The first and second derivatives of u\(£,t) enter the problem also. The corresponding equations can be obtained by 
taking derivatives of Eq. (|V.ll[) for s = 1. 

The recursion is consistent: Jacobi fields <9*i(*(£, t)/c?£ !i ...<9£ lt are calculated sequentially towards increasing orders 
s and t. The derivatives of Eq. (|V.ll[) involve derivatives of u\ (£, r) with r < s only. 

In the potential scattering of four particles (n = 12), we would have to calculate 24, 7800, and 499200 functions 
describing characteristics and Jacobi fields to orders h°, h 2 , and ft 4 , respectively. As compared to the previous method, 
we gain a reduction in computer time by several orders of magnitude. In the potential scattering of two oxygen atoms 
le O + le O (16 electrons and two nuclei, n — 54 degrees of freedom), one has to calculate 108, 647460, and ~ 10 10 
functions of time describing characteristics and Jacobi fields to orders h°, h 2 , and h 4 , respectively. The second order 
calculation might numerically be feasible. 
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C. Monte Carlo method for averaging over the Wigner function 



The average value of an observable associated to an operator f can be calculated using equation /(£, r) = 
r), 0), expansion (|IV.7p . and the initial-state Wigner function W(t;,0) = W(t;): 

f d 2n £ 

< m, r) >= J j—fi;f(*u(Z, r))W(0. (VII.3) 

We describe calculation of (IVII.3|) using an extension of the Monte Carlo method exploited in lattice QCD [65]. 

Let us separate phase space into two regions f2+ and £!_ where the Wigner function is positive and negative definite, 
respectively. One has W(£) = — where VF±(£) > 0. The next step deals with generating events in f2± 

distributed according to probability densities W±(£). 

Consider region f2+. One starts from generating 2n + 1 numbers (£ J ,7) with £ l distributed homogeneously in Q + 
and 7 distributed homogeneously in the interval (0, W max ) where 

W mm =max.W+{0. (VH.4) 

The joint probability density of the 2n + 1-dimensional variable (£ J , 7) is assumed to be proportional to 8(W+(£) — 7). 
One can check that, owing to a constant normalization factor, the marginal probability density of £ l is given by 

W+(0= 9(W+(0-l)dj. (VII.5) 



Given (£',7), we examine condition 



W+(£)-7>0. (VII.6) 



If it is fulfilled, we shift the number 2V+ of successful events by one unit, set a — £, calculate quantum characteristic 
it l (£+ a ,r), observable f(*u(£ t+ a , r)), and keep the record with a = N + . The star-product from the arguments of 
/(*«(£+ q ,t)) can be removed as described in Sect. IV. If inequality (|VII.6[) is not fulfilled, next 2n + 1 numbers 
(£ S j7) are generated. 

The similar procedure applies for f2„. 

Suppose we have generated samples of the iV+ and iV_ successful events £± a G fi± ■ To get the average value of 
function /(£, r), the values f(*u(£± a > T )) should be multiplied by 

w ± = I -Br wm - (VIL7) 

where — W- — 1, divided to total numbers N± of successful events, and summed up to give 

< m,r) >« ^E/MC+^r)) -^Y^f(*u(£_ a , T )). (VII.8) 

+ a=l a=l 

This equation accomplishes reduction of the evolution problem of an n-body quantum-mechanical system to 
a statistical-mechanical problem of computing an ensemble of quantum characteristics u l (^± a ,r) and gradients 

d t U i (S± a ,T)/d?K..d?*. 

The Monte-Carlo method is efficient for calculation of multidimensional integrals of real positive functions. In 
quantum mechanics, the amplitudes are complex functions which oscillate rapidly in the semiclassical regime. The 
numerical calculation of oscillating functions requires the high precision, diminishing therefore the value of numerical 
methods for calculation of path integrals in the configuration space. The imaginary time formalism combined with the 
Monte-Carlo method is effective for finding binding energies of the lowest states, but not for scattering problems. The 
path-integral method in the phase space (oO . fsij ] represents perhaps more promising tool, since the Green function is a 
real function, it is positive definite in the classical limit, although it is not positive definite in general. The phase-space 
path integral method is used in Ref. [66] to find the evolution of a quantum state as an illustrative example (see also 
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D. Scattering theory in terms the Wigner function 

The scattering theory treats elementary and bound particles essentially on the same footing. We fix in an out 
scattering states at r = t' and r = t" , respectively, and select clusters a and (3 for elementary and bound particles 
with definite momenta and 

W m {0 = Y[^(2irhfS(p' a -Y.P^ W ^> ( VIL9 ) 

Here, £ = (xi, xjv, Pi, Pat), N is the total number of particles participating in the reaction, the number of 
clusters is less or equal to N, V is the normalization volume. The Wigner functions Wj n (£) and W ou t(C) correspond 
to products of free plane waves of elementary particles and bound states and include the product of the intrinsic 
Wigner functions of bound states. We thus have W = W Q (£ Q ) = 1 for elementary particles, whereas for bound 
states Wp(£p) and W a (£ a ) represent the Wigner functions of bound states in the rest frame, which are supposed to 
be known. The transition probability from the initial in-state in to the final cwi-state is the absolute square of the 
S-matrix element Wfi = | < out\in > | 2 . In terms of the Wigner function, one finds 

f d 2n £ 

*»fi = J j^w out (M^r))w in (0, (vn.ii) 

where n = 3N. At this stage, Eq.(7.8) can be used for numerical simulations. 

Equation (jVII.lip simplifies for scattering of one particle in an external potential. The cross section is defined by 

da ji = - — lim lim ' 



where ji nc — tf- In the classical limit 



jinc t"^+mt'^-oo t" — T 1 (2ttH) 3 



± J dx(2 7 rft) 3 ( 5 3 (p"^^+ l (x,p',r)) 



The initial momentum p' is directed along the z-axis. The scattering takes place on a finite-range potential. The limit 
for the integral along the particle trajectory extends from z = to z — v'(t" — r') — > +oo owing to contributions O(l) 
coming from the interaction region, which can be neglected. The integrand does not depend on z, so this contribution 
is factorized to give simply v'(t" — r'). As a second step, the change of variables (0, b) — > (</>, 9) can be done, where <j> 
is the azimuthal angle, b is the impact parameter and 9 is the scattering polar angle. In this way, one gets the well 
known expression 

b db 
da c f [ = — dOp«. 

The semiclassical expansion of W ou t(*u(^, r)) around W ou t(u(t;, r)) can be useful to calculate quantum corrections 
to the classical cross sections. 



VIII. CONCLUSIONS 



In this work, we discussed properties of the Weyl's symbols of the Hcisenberg operators of canonical coordinates and 
momenta. Functions t) which show quantum phase flow obey Hamilton's equations in quantum form (|III.11[) . We 
derived these equations using the Green function method. The principle of stationary action for quantum Hamilton's 
equations was formulated and quantum-mechanical extensions of the Liouville theorem on conservation of the phase- 
space volume and the Poincare theorem on conservation of 2p-forms were found. 

Symbols of generic Heisenberg operators are entirely determined by r) from 



= /(*«(£, r),0). 
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This equation is remarkable in many respects. It shows that quantum phase flow comprises the entire information 
on quantum evolution of symbols of Heisenberg operators. Usually characteristics satisfy first-order ODE e.g. the 
classical Hamilton's equations and solve first-order PDE e.g. the classical Liouville equation. Functions w l (£, r) are the 
genuine characteristics of Eq. (|III.14[) . despite both m j (£,t) and /(£, t) obey infinite-order PDE. It is quite surprising 
that to any fixed order of the semiclassical expansion w j (£,t) and therefore /(£, r) can be constructed by solving a 
system of first-order ODE only. Such a circumstance allows to approach the problem numerically using efficient ODE 
integrators. 

We discussed methods of eliminating the ^-symbol from arguments of functions /(*u(£,r),0) by semiclassical 
expansion around /(u(£, r), 0) and provided the explicit expressions up to the fourth order in h. 

The energy conservation in quantum form 9[) renders upon expansion over the Planck's constant an infinite 
sequence of conserved quantities Eq. (|V.14[) depending on u*(£,t) and finite-order derivatives of u 1 s (£_,t) with respect 
to i.e. Jacobi fields. The order h 2 conserved quantity (|V.18[) is expressed in terms of classical trajectories and their 
Jacobi fields only. 

For illustration purposes, we constructed the lowest order quantum corrections to the Kepler periodic circular orbits. 
These corrections increase linearly with time and show the resonance behavior. 

The phase-space formulation of quantum dynamics can be useful for transport models in quantum chemistry 
and heavy-ion collisions where embedding coherence in propagation of particles remains, in particular, an unsolved 
problem. The principal advantage of the deformation quantization is its proximity to the classical picture of phase- 
space dynamics and the strict quantum character concurrently. Specific effects such as quantum coherence and 
non-locality are displayed implicitly upon the H expansion in the dependence of physical quantities on Jacobi fields. 

The main result of this work is the reduction by semiclassical expansion of quantum-mechanical evolution problem 
to a statistical-mechanical problem of computing an ensemble of quantum characteristics and their Jacobi fields. 
The reduction to a system of ODE pertains rigorously at any fixed order in h. Given quantum characteristics are 
constructed, physical observables can be found without further addressing to dynamics. The method of quantum 
characteristics can be useful for calculation of potential scattering of complex quantum systems. 

This work is supported in part by DFG grant No. 436 RUS 113/721/0-2, RFBR grant No. 06-02-04004, and 
European Graduiertenkolleg GR683. 
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